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Abstract 

We deal with a boundary detection problem arising in nondestructive testing 
of materials. The problem consists in recovering an unknown portion of the 
boundary, where a Robin condition is satisfied, with the use of a Cauchy 
data pair collected on the accessible part of the boundary. We combine a 
linearization argument with a Tikhonov regularization approach for the local 
reconstruction of the unknown defect. Moreover, we discuss the regulariza- 
tion parameter choice by means of the so called balancing principle and we 
present some numerical tests that show the efficiency of our method. 

Keywords: Tikhonov regularization, Robin boundary condition. Free 
boundary problem. Balancing Principle, Local identification 



1. Introduction 

In this paper we deal with an inverse problem arising in corrosion detection. 
We consider a domain C which models a 2D transverse section of a thin 
metallic specimen whose boundary is partly accessible and stays in contact 
with an aggressive environment. Hence, in order to detect the damage which 
is expected to occur in such a portion of the boundary, one has to rely on 
the electrostatic measurements of a potential u performed on the accessible 
portion. 
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We are then lead to the study of the following elliptic boundary value problem 



Am = , in f2 

— = $ , on 

— + 7M = , on T/ 
ov 

M = , on r^i 



According to this model u is the harmonic potential in VL. We assume that the 
boundary of Vt is decomposed in three open and disjoint subsets Va, T/, Vr,. 
On the portion F^, which is the one accessible to direct inspection, we pre- 
scribe a current density $ and we measure the corresponding voltage poten- 
tial u\yj^- The portion F/, where the corrosion took place, is out of reach. 
On such a portion the potential u satisfies an homogeneous Robin condition, 
which models a resistive coupling with the exterior environment by means of 
the impedance coefficient 7. 

In this paper we are interested in the numerical reconstruction issue of the 
unknown and damaged boundary F/ from the data collected on the accessible 
part of the boundary F^, that is the Cauchy data pair {u\t^,^). 
Boundary and parameter identification results related to this stationary in- 
verse problem has been provided by many authors [H El El El |9l IH [THl ttH fT2| 
[131 [a [151 HSl [IT]. 

Local uniqueness and conditional stability results for the inverse problem at 
hand are contained in ^ and constitute the theoretical setting on which our 
numerical analysis relies. The present local determination of corroded bound- 
aries consists in the localization of a small perturbation Tj^q of a reference 
boundary F/. It is convenient to introduce a small vector field 9 G Co(F/) so 
that the damaged domain VLq is such that 



dVte = Fa U Fi) U F/^e, Vj^e = {z : z = w + 9{w), w G Tj}. 

Such a local approach combined with a linearization argument (see [5]) allows 
a reformulation of the problem of the localization of the unknown defect Tj g 
as the identification of the unknown term 6' in a boundary condition of the 

type 

du' , d d \ ^ , 

+ ^u' = — [9- v—u + 7^ ■ z/ 7 + 2H)u 
Of as \ as 
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at the portion Tj, where u' is a harmonic function satisfying homogeneous 
Neumann and Dirichlet conditions on Ta and Td respectively, u is the solu- 



tion of (1.1), and H denotes the curvature of the reference boundary Fj. As 
in ^ we carry over our analysis under the a-priori assumption of a constant 
7 such that 2H{x) + 7 > in F/ and we limit ourselves to the case of positive 
fluxes $ only. 

We linearize the forward map F : ^ H- Melr^) where ug is the solution of the 



system (2.3) below, by its Frechet derivative F' and take the voltage contrast 



on Ta, as the noisy right-hand term for the considered operator equation, 

F'e = {ue-u)\r^. 

As in ^13j, we assume that the unavoidable measurement errors in voltage 
contrast are not smaller than the truncation error, o(||^||). Therefore, if the 
noise level for voltage measurements is assumed to be 6, then the noise level 
for the right-hand term of the above operator equation can be written as 
6 = K6, where a constant K is not necessary to be precisely known. Our 
method is based on a discretized Tikhonov regularization argument where 
the regularization parameter is chosen by a balancing principle (cf.[6| [8| 120]). 
Such an a posteriori parameter choice can lead to a regularized solution with 
order-optimal accuracy. At the same time it can provide a reliable estimate 
for the constant K. 



2. Local identification of the unknown boundary 

In this section we shall collect the main identifiability results of which our 
reconstruction procedure and our numerical tests are a follow up. For a more 
detailed description we refer to [Sj. 

We denote with u the outward normal to F/ and we assume that 6' is a vector 
field in C^{ri) having a nontrivial normal component 6^ on F/. 
Let the Sobolev space Hq{Q,T£)) be defined as follows 

H^{n, Tn) = {v e H^{n) : v = on in the trace sense} . (2.1) 

We introduce the forward map F 

F: C',{Tj) ^ H^Ta) (2.2) 
6 ^ ue\rA 
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where ug G Hq{Q,Td) is the solution to the elhptic problem 

in Qg 

on Ta 



Aug 
dug 



OUg ^ 

du 

u = 



■ $ 



on Tj g 
on Yd- 



(2.3) 



We recall the following differentiability property for the forward map F. 



Lemma 2.1. The operator F in (2.2) is Frechet differentiable at zero. In- 
deed, consider the linear operator F' : Co(r/) — H^ITa) defined as F'6 = 
u'\r^, where u' is the solution to the boundary value problem 



( Au' 

du 
du' 



u 



du 
' - 



+ 7^' 



d 
d^ 



in Q 
on Ta 



e^—u ) +70^(7 + 2i/)M on Tj 



(2.4) 



on r 



D, 



the function u is the solution of (1.1) and H denotes the curvature of the 
boundary Tj. Then, 



\\F{e) - F(0) - F'eW^r^^^^ ^0 ase^0^n C^iTj). 



Let us also recall that a weak solution to (2.4) is a function u' G HqIQjTd) 
such that 



Vw'Vf + 



'JU V 



76^^(7 + H)uv 



d d 
as as 



(2.5) 



for all V G f^o^(^], Fz?). 

The following theorem ensures that the operator F' is injective, under some 
reasonable hypothesis. This property allows us to conclude that the solution 
6 to our inverse problem is identifiable, at least for small perturbations. 
Moreover, we recall a conditional Lipschitz type upper bound for on a 
suitable portion of Tj in terms of m'Ita = thus showing that the inversion 
of F' is not too much ill-behaved. 
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Theorem 2.2. Let $ e ^^5(1^) be nonnegative in the sense of distributions. 
Let us assume that 2H{x) + 7 > and 6y{x) ^ for any a; G F/. Then F' is 
injective. Moreover, there exists a positive constant c > such that 

\\u'\\ 1, f \9\ 

II ii/f2(rA) ^'11 



where Tj is an inner portion of the boundary Tj. 

Finally, in the next theorem, we consider L^(Fyi) as codomain space of the 



operator F' introduced in Lemma 2.1 stating a compactness result. 



Theorem 2.3. The linear operator 



F': Cl{T,) L\Ta) 



(-7- u\y 



A 



where v! is the solution to the boundary value problem (2.4), is compact. 



3. Tikhonov regularization for a local reconstruction and an esti- 
mate of the accuracy 

Here and in the following, with a slight abuse of notation, we shall denote 
by F' the compact operator 

F' : El{V,) ^ L^(F^) (3.1) 



where u' satisfies the weak formulation in (2.5) for any v G ifo(fi,F£)). 
The existence and uniqueness of v! G H\(VL,Vd) follows from standard ar- 
guments on elliptic boundary value problems. Moreover, the compactness of 



F' in (3.1) follows along the lines of Theorem 4.5 in 

In view of this compactness property, the issue of the identification of Q may 
be interpreted as the regularized inversion of the above compact operator 
F' = F'(Tj) between the Hilbert spaces Hq(Ti) and L^(Fa). Such kind 
of reformulation allows us to deal with the approximate inversion by the 
technique of Tikhonov regularization. 

We are interested in finding the solution to operator equation 

F'e = f:=u'\r^, (3.2) 
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where instead of the exact data r, a noisy version is known. As in 



if we hnearize the forward map F defined in (2.2) by its Frechet derivative, 



then by Lemma 2.1, we obtain 



I.e. 



F'e = F{e)-F{o) + om\), 



F'e = iue-u)\r^+om 



Here ug\r^ and -ulr^ are voltages measured in experiments. In practice they 
are usually given in a noisy form as ug\^^ and m|p^ with 6 being the noise 
level for unavoidable experimental error for the measurements of the voltage. 
When II ^11 is rather small, one can assume that these measurement errors 
in voltage contrast [uq — u)\rj^ have the same order of magnitude as the 
truncation error o(||6'||). Thus, we take 

/ ■.= ul\r,-u'\r. 



as the noisy right hand term for (3.2) and assume that 

C6, 



< 6 



(3.3) 



where a constant C is unknown. 

If the Tikhonov regularization is applied to the ill-posed operator equation 

F'e = r\ 
then the regularized solution solves 

{F')*F'e + al= {F')*r^ (3.4) 
where a > is the Tikhonov regularization parameter and / is the identity 



operator on space Hq^Ti). It is well known that the solution to (3.4) will be 
the minimizer of the functional 



J{e) := ||F'^ - r'WUrA) + «ll^llHi(r,)- (3-5) 
Here, we assume the exact solution 9 belongs to the set of source condition 
Mh ■■= {s e H^iTj) : s = h{{F'yF')w, \\w\\ < l} (3.6) 
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where h is an 'index function' defined on [0, oo), which is operator monotone 
(see [iHl [19]) and satisfies the condition h{0) = 0. Moreover, it has been 
proven that 



sup 

0<A<6 



a 



a + \ 



h{X) < h{a) for all a e (0, a] (3.7) 



and some a > 0. 

Let us notice that J^{0) in (3.5) is the standard Tikhonov regularization 
functional, where the penalty term naturally is imposed in i^Q-norm. Such 
a consideration can facilitate the analysis for the accuracy. Moreover, it is 
equivalent to the Tikhonov regularization functional considered in [13] with 
a penalty term based on the derivative of the regularized solution. 



The discretization of the regularized problem (3.4) is realized by Galerkin 
method. The Galerkin approximation of Tikhonov-regularization consists 
in minimizing the above functional J^{0) in a finite-dimensional subspace 
Xn C Hq(Tj). As usual, in Galerkin scheme, the discretized regularized 
solution 9^, „ is characterized by the variational equations 

{F'el^ - r\ F'z) + a{el^, z) = 0, G (3.8) 

or, equivalent ly, 

where := -F'P„, with P„ being the projection from HliVi) onto X„. 
Let /i, /2, . . ., /„ be basis functions of X„. If one decomposes 6^^^ into a 
linear combination of /i, /2, . . ., /n, i-e. 0^^ = Yl^=i ^ifii then the coefficient 
vector c = {q}"^]^ can be obtained by solving a linear algebraic system, 

{M + aG)c = Rs, 

with the following matrices and vector, 

G := ((/.,/,)H„Hr,))". ^ (3.10) 
Rs := {(F7.,r^)L^(r,)}L,. 



Remark 3.1. The adjoint operator {F')* is not involved in the construction 

a,n 



of 9^„ . Theoretically, F'fi can be obtained by solving the boundary value 
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system (2.4) and deriving the trace on Ta, where function 9 is replaced by 
fi, for i = 1, . . . ,n. Moreover, we do not need each F'fi in an explicit form, 
but only its products in (3.10), which can be computed much more accurately 
than F'fi itself. 

According to the classical results on Tikhonov regularization for linear ill- 
posed problem and in view of (3.3) and (3.7), it holds that 

1 1 ZD nS II 

II a.«ll/fi(rj) 



< 



+ 



d-nKrK+o^iy'iKrr ^^^^^^^ 



< 



e-{{F'^rF'^+air\F'^rf 



C5 

+ 



As in [19], we can estimate the noise free term as follows, 
e-{{F'J*F'^ + aI)-\F'J^f 



< 



< 



+ 

|( 
+ 

+ 



{{F'J^F'^ + aiy\F'J*F'{I-P^)e 



[i - iiKTK + ociy' iKTK) ihiiF'TF') - h{{F':)*F':)) w 



a 



< Ci{h{a) + h(\\F'{I -P, 



\F' {I-Pn 



where the constant Ci does not depend on a and n. 

In view of the best possible order of accuracy without discretization being 
h{a) + 5/ a/o, the discretization has to be chosen such that 

(/ - Pn) : Hl{Ti) ^ L\Ta)\\ < 6. (3.11) 

Summing up the estimates above, we have the following theorem. 
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Theorem 3.1. Under assumptions (3.3) and (3.6), and with discretization 

5 



satisfying (3.11), there holds that 



\9-9iJ <Kh{a) + K- 



(3.12) 



where the constants K and K do not depend on a and 5. 

4. Parameter choice rule based on the balancing principle 

In this section, we give a regularization parameter choice rule based on the 
balancing principle developed in [6l [71 [20]. The essential idea of this prin- 
ciple is to choose the regularization parameter a balancing the two parts in 



error estimate (3.12). As a posteriori parameter choice rule, the balancing 



principle can select regularization parameter in an adaptive way without a 



priori knowledge of the solution set (3.6). That is, the index function h in the 



bound (3.12), which indicates the smoothness of 9 as shown in (3.6), does not 



need to be known. At the same time, it does not require to know the precise 



noise level, either. In our model problem, constant C in (3.3) indicating the 



precise noise level is unknown, which leads to K in (3.12) is also unknown. 



A reference noise level 5 is sufficient for the performance of the balancing 
principle. The regularization parameter chosen by the balancing principle 
leads to a regularized solution with an order-optimal accuracy. 
Assume that the projection P„ is chosen with n 
is satisfied. Let 9^^ := 9l^^^^^p,^^y 
We select parameter a from the geometric sequence 



n{F\5) such that (3.11) 



A:={a„ = aog", n = 0, 1, . . . , A^}, 
with g > 1, sufficiently small oq, and sufficiently large N such that aN-i < 

1 < ttAT. 

For any given K , one can choose the parameter from A(ao,g, A^) by the 
following adaptive strategy. 



a{K) = max 



a„ G A : 11^^ - 9i ll^ifr ^ < K6 + 



m = 0, 1, 



,n — 1 



We further rely on the assumption that a two-sided estimate 



(4.1) 



(4.2) 
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holds for some c G (0, 1), where 6^^ is defined by (3.9) with being taken as 
f. The upper estimate in (4.2) is due to (3.12). As to the lower estimate, 
it just suggests that the noise propagation error is not that small. If the 
lower estimate is not satisfied, it just means that our estimate to noise level 
is too pessimistic. However this will not cause a problem, since later we 



shall show that under assumption (4.2) balancing principle can provide an 
order-optimal accuracy. 

Now, consider the following hypothesis set of possible values of the constant 
K 

lC = {kj = koP^, j = 0,l,...,M}, p>l, 
and assume that there are two adjacent terms ki, fc^+i G /C such that 



ki<cK <K < k+i. 



(4.3) 



In fact, each element in /C can be viewed as a candidate for the estimator to 
constant K. Our aim is to detect fc^+i (or say ki) among the elements in JC, 



and to use ki+i in adaptive strategy (4.1) to obtain a parameter a. 

In view of (4.2) and (4.3), if the hypothesis kj G JC for K is chosen too small. 



i.e., kj < ki then, as it is shown in [6l [20], the corresponding regularization 
parameter a{kj) will be smaller than a threshold depending on and p. 
Thus, if 



a{ki) := min < a{kj) > 9ao 



p- 



(4.4) 



then there holds that, either z = /, or i = / + 1. 

In order to guarantee the regularized solution is stable enough, we choose 
final regularization parameter as 

a+ = a{ki+i). 

With such a choice a+, we have the following theorem. 
Theorem 4.1. Under the assumptions above, there holds 

\\0 - OijHi^r,) < Qp'VQKhCh~\K6)), 



where h{t) = Kh{t)\/t, K is the constant from estimation (3.12), and h ^ is 
the inverse function of h. 
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Note that from 



it follows that the error bound indicated in Theorem 4.1 

2 



is order-optimal, i.e., it is only worse by a constant factor 3p\/q than a priori 
optimal bound 2Kh{h~^{KS)). If index function h in the source condition 
(Q is given as h{X) = cX", < u < 1, then h{h-\K6)) = 0{6^), which 
coincides with the classical rate for Tikhonov regularization. 



Remark 4.1. The proof of Theorem 4-1 can be referred to or For 
the general discussions on the application of the balancing principle with two 
flexible parameters, one can refer to J2l 



5. Numerical tests 

In this section, we present some numerical examples to illustrate the theo- 
retical results obtained above. 



In Examples [IH3| we consider the corrosion problem in (1.1) with 

Vl = (0,7r) X (0,1), 

rA = (0,7r)x{0}, 

r, = (0,7r)x{l}, 

r^ = {0}x (0,l)U{7r}x (0,1). 

On such a rectangle domain, if the flux $ = sin(a;) is given on F^, then the 



solution to (1.1 ) is in the form of 



/ X f • , / X 7sinh(l) + cosh(l) , \ . / n 

u{x,y) = -smh(y) + . ——cosh?/ sm(a;) 

\ smh(l) + 7cosh(l) / 

with 7 > 0. We test the same flux $ = sin(a;) in Examples [l]|3} 

Example 1. 

In this example the vector fleld 6 is given as ^ = (0, ^2(2;)) with 



eoix) 



e'2{t)dt, where ^2(0) = 0, 



and 



9'ix) 



cot(x) + 7W cot^(x) —7^ + 1 



cot^(a;) 



7 



TT 

-, <x < -, 



— cot(x) — 7-v/cot (x) —7^ + 1 



cot^(a;) — 7 
11 



7r 



with the constant 7 such that < 7 < 1. Such a choice of 9 corresponds to 

ue{x,y) = exp{-y) sin(x) 



solving (2.3) 




(a) 02 (x) 



(b) 9',ix) 




1,5 2 2,5 3 




(c) noise added to {uq — u)|rA (d) comparison of r = u'|r^ and 
Figure 1: Functions in Example [T] with 7 = 0.999. 

As we mentioned in the Introduction, the impedance coefficient 7 depending 
on the exterior environment should be a fixed constant in the model problem 
(1.1). However, in this particular example, the scale of 6 depends on 7. 



On the other hand, our linearization approach by truncation can only work 
when ll^ll is rather small. Thus, in this example we test different values of 7 
which are all quite close to 1. Figures [l] (a) and (b) illustrate the behaviors 
of ^2(2;) and 02{x) when 7 = 0.999. In order to simulate the error arising 
in experiment measurements, we add random noise to each grid involved in 
calculation, i.e. we take 

where ^ is random variable with range [—1,1] and the reference noise level 6 = 
10~^. Figure [l](c) and (d) show the additional noise ^6 and the comparison 
of f = u'Ita ^iid in the case 7 = 0.999. 



In all of the following tests, the discretization level in (3.8) is taken as n 



20 and the regularization parameter a is chosen by the balancing principle 
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described in Section |4j In the case of 7 = 0.999, the parameters in the 
implementation of the balancing principle are settled as follows: 



A 



-- aoq , n 
koP^, 3 = 



= 0,1, 
0,1,.. 



■,N}, 
M}, 



ko -- 



--1-10 
0.006 



-11 



P 



q = 1.3, N = 69] 
= 1.3, M = 19. 



Sequence {a{kj)}^^Q produced by (4.1) with K replaced by kj results in 



4.827- 10-", 
9.173- 10-^ 
1.644- 10^^ 
2.946 - 10-^ 



8.157- 10-", 
3.406 - 10-^ 
2.778 - 10-^ 
8.415 - 10-^ 



1.379 - 10-1°, 
7.482 - 10-^ 
3.612 - 10-^ 
1.094 - 10-5, 



3.937-10-1°, 
9.728 - 10-^ 
6.104- 10-^ 
1.422 - 10-5, 



1.900- 10-^ 
1.265- 10-^ 
1.341 - 10-^ 
1.849 - 10-5. 



For the parameters designed as above, the value of the threshold is calculated 
2 



as 9ao 



p-i 



7.236 - 10-^ Then a. 



time, we obtain an estimate to K as ky 



aikj) = 3.406 - 10-^. At the same 
= 0.0290, which suggests the true 



noise level KS is about 2.90 - 10-*^. 

Table [T] summarizes the results for the other values of 7 in Example [T] 
Here we take different reference noise level 6 according to 7, because, as 
we mentioned, in this particular example, 7 determines the scale of 6 and 
furthermore the truncation error. In Table [l| Ervhi '■= ||^ — da^^m{Ti) and 
11^^ — 6a^\\L'^{Ti) denote the errors in the corresponding norms. The 



Err, 
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reconstructed functions 9i , are displayed in Figure |2 



Example 2. 



In this example, the vector field 6 = {0,62{x)) to be identified is similar 
to what is considered in [T3], where ^2(2^) is a piecewise linear function, 
as shown in Figure |3| In contrast to Example [T| we do not assume that 
7 G (0,1), and test two cases: 7 = 1, 7 = 10. The solution of (2.3) and 



7 


6 


K 


a. 




Errhi 


Erri2 


0.9999 


io-« 


0.0489 


2.015 - 


io-« 


0.0012 


9.134 - 10-5 


0.999 


10-6 


0.0290 


3.406 - 


io-« 


0.0080 


6.127- 10-^ 


0.99 


10-5 


0.0636 


1.265 - 


10-^ 


0.0438 


0.0038 


0.95 


10-^ 


0.0489 


5.756 - 


10-8 


0.1482 


0.0308 



Table 1: Test results for Example [T] 
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(a) 7 = 0.9999 



(b) 7 = 0.999 



a=1.265 10"' 



cc=5. 756-10"' 





(c) 7 = 0.99 



(d) 7 = 0.95 



Figure 2: The simulated solutions 9^^ in Example [l| 



: =5.756-10"' 



: =9.728 10"' 





0.5 1 1.5 2 2.5 3 



(a) 7 = 1 (b) 7 = 10 

Figure 3: The simulated solutions 6^^ in Example [2| 



7 


6 


K 








Erri2 


1 


10-7 


0.3937 


5.756 ■ 


10-8 


0.0031 


3.565 ■ 10"^ 


10 


10"7 


0.5119 


9.728 ■ 


10"8 


0.0030 


2.622 ■ 10"^ 



Table 2: Test results for Example [2] 



0.5 1 1.5 2 2.5 3 



0.5 1 1.5 2 2.5 3 



its trace uelr^ ^I's generated numerically. We simulate point-wise random 
noise in each discretization note on with reference level S = 10"^, and 
a+ is chosen according to the balancing principle under such a value of 6. 
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The approximations 6^^^ are displayed in Figure [s] and the test results are 
summarized in Table [21 

Example 3. 

In this example, we take the vector field 9 = {0,62{x)) with 



e2{x) = h- 




< X <TT, 



as shown in Figure |4j Here one can change the value of the constant /i > to 




(a) h = SOtt (b) h = 907r 



Figure 5: The simulated solutions 6^^ in Example [sj 

adjust the scale of 6. The solution U0 and its trace on in this example are 
also obtained numerically. In order to guarantee that the scale of 6 is small 
enough such that the truncation method can work well, we test h = SOvr and 
h = QOtt. In both cases, larger value of 7 may make the problem less ill-posed 
and result in better reconstruction. The approximations 6'^^ are displayed 
in Figure [5] and the test results are summarized in Table |3} 
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h 


7 


6 




K 


«_ 






Ervhi 


Erri2 


SOtt 


1 


10" 


-5 


0.0371 


1.341 ■ 


10" 


-6 


0.0079 


0.0034 


10 


10" 


-5 


0.0371 


1.849 ■ 


10" 


-5 


0.0016 


4.273 ■ 10-4 


gOvr 


1 


10" 


'7 


0.1379 


3.611 ■ 


10" 


-6 


0.0027 


0.0013 


10 


10" 


-7 


0.1060 


1.743 ■ 


10- 


-6 


6.300 ■ 10-4 


1.340 ■ 10-4 



Table 3: Test results for Example [3] 





3 




^0 







Figure 6: Domain 0.0 in Example |4] 



Example 4. 



In the last example, we consider a domain Vt given as a half annulus bounded 
by the following curves (see Figure [6]). 



(x, y) : y = V4 - x2, -2 < x < 2 , , 

{x,y) : y = Vl - a;2, -1 < x < l| , 
rz) = {(x,y): y = 0, l<|x|<2}. 



For flux $ = y on F^, the function m solving (1.1) can be written as 



7 


5 


K 


a. 


\- 


Ervhi 


Erri2 


0.99 


10-5 


0.0816 


1.032- 


10-5 


0.0057 


0.0028 


0.9 


10-4 


0.0816 


1.032- 


10-5 


0.0708 


0.0384 



u{x,y) =Ay + 



Table 4: Test results for Example [i] 
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exact solJtion 
simulated solution 




(a) 7 = 0.99 (b) 7 = 0.9 

Figure 7: The simulated solutions in Example [ij 



with < 7 < 1. For vector field 9 = {0,62{x)), we consider 62{x) 
\/l — x^, where 



if{x) 



^{x)={ J 



-a/1 - (7X + 7- 1)^ 



-1 < X < 0, 



7 



l-(7X-7 + l)2, 0<x<l. 



It can be verified that ug{x,y) = y solves(2.3) in ilg. 

In this example 7 also determines the scale of 6. Thus, we take the values of 
7 very close to 1. The approximations 9^^ are displayed in Figure [t] and the 
test results are summarized in Table HI 

We would like to note that in all considered examples the balancing prin- 



ciple (4.1), (4.4) has been implemented with the same values of the design 



parameters ao, p and g, because the domain f2 and the operator F' are the 
same for all examples. This suggests that in practice, for a given domain 
Q the parameters Oq, p and q can be determined in the experiments with a 



problem (2.4) where a solution is known, and then kept for studying all other 



problems (2.3) in the given domain f2. 



Both the theoretical and numerical results suggest that the linearizaton ap- 
proach considered in this paper can perform well for the identification of the 
corroded boundary only on condition that the scale of this boundary function 
is quite small. This is the limitation of the approach. However, in practice 
one certainly does not expect too much corrosion taking place to the metallic 
specimen. 
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